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We have carried out a wide study of Cold Nuclear Matter (CNM) effects on J/tf/ production in dAu , CuCu 
and AuAu collisions at y%jv = 200 GeV. We have studied the effects of three different gluon-shadowing 
parametrisations, using the usual simplified kinematics for which the momentum of the gluon recoiling against 
the J/if/ is neglected as well as an exact kinematics for a 2 —> 2 process, namely g + ^ — > J/iff + g as expected 
from LO pQCD. We have shown that the rapidity distribution of the nuclear modification factor R^Am and 
particularly its anti-shadowing peak, is systematically shifted toward larger rapidities in the 2 — > 2 kinematics, 
irrespective of which shadowing parametrisation is used. In turn, we have noted differences in the effective 
final-state nuclear absorption needed to fit the PHENIX dAu data. Taking advantage of our implementation of 
a 2 ^ 2 kinematics, we have also computed the transverse momentum dependence of the nuclear modification 
factor, which cannot be predicted with the usual simplified kinematics. All the corresponding observables have 
been computed for CuCu and AuAu collisions and compared to the PHENIX and STAR data. Finally, we 
have extracted the effective nuclear absorption from the recent measurements of Rcp in rfAu collisions by the 
PHENIX collaboration. 

PACS numbers: 13.85.Ni,14.40.Pq,21.65.Jk,25.75.Dw 



I. INTRODUCTION 

The J/i// particle is considered to be one of the most in- 
teresting probes of the transition from hadronic matter to a 
deconfined state of QCD matter, the so-called Quark-Gluon 
Plasma (QGP). In the presence of a QGP, binding of cc pairs 
into J /if/ mesons is predicted to be hindered due to color 
screening, leading to a suppression in heavy ion colli- 
sions [1]. 

The results of the J/tfr production in AuAu collisions at 
-\fsNN-2QQ GeV [2] show a significant suppression. Never- 
theless, PHENIX data on dAu collisions [3] have also shown 
a non trivial behaviour, pointing out that Cold Nuclear Mat- 
ter (CNM) effects play an essential role at these energies (for 
recent reviews see [4—6] and for perspectives for the LHC 
see [7]). 

All this reveals that the interpretation of the results obtained 
in nucleus-nucleus collisions relies on a good understanding 
and a proper subtraction of the CNM effects, known to impact 
the J/ip production in deuteron-nucleus collisions, where the 
deconfinement cannot be reached. 

In previous studies [8, 9], we have shown that the impact 
of the gluon shadowing on J /if/ production does depend on 
the partonic process producing the cc and then the J /iff. In- 
deed, the evaluation of the shadowing corrections in which 
one treats the production as a 2— »2 process {g + g ^ J/tf/ + g) 
shows visible differences in the nuclear modifications factors 
when compared to the results obtained for a 2— >1 process. 

The former partonic production mechanism seems to be 
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favoured by the recent studies based on Colour-Singlet Model 
(CSM), including [14] or not [15] i-channel cut contribu- 
tions. This also seems to be confirmed by the recent stu- 
dies of higher-order QCD corrections, which have shown, 
on one hand, that the problematic Pj dependence of the 
LO CSM [17] is cured when going at or^ and [18-22] 
and, on the other, that the CSM yield at NLO in e+e" -> 
J/i// + Xnon cc [23, 24] saturates the experimental values by 
the Belle collaboration [25]. This does not allow for a sig- 
nificant Colour-Octet (CO) component [26], which happens 
to be precisely the one appearing in the low-f 7 description of 
hadroproduction via a 2 — > 1 process [27]. To summarise, one 
is entitled to consider that the former 2 — > 2 kinematics is the 
most appropriate to account for the PHENIX data. 

The structure of this paper is as follows. In section 11, 
we describe our approach, namely the partonic process, the 
shadowing parametrisations and the implementation of the 
nuclear absorption that we have chosen. In section III, we 
present and discuss our results for Rjau versus rapidity, cen- 
trality and transverse momentum. We particularly discuss the 
impact of the partonic process kinematics. Section IV is de- 
voted to the results in AuAu and CuCu collisions. In section 
V, we present and discuss our extraction of the effective J/iJ/ 
absorption cross section from the PHENIX dAu data and we 
finally conclude. 



II. OUR APPROACH 

To describe the J/i// production in nucleus collisions, our 
Monte-Carlo framework [8, 9] is based on the probabilistic 
Glauber model, the nuclear density profiles being defined with 
the Woods-Saxon parametrisation for any nucleus A > 2 and 
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the Hulthen wavefanction for the deuteron [12]. The nucleon- 
nucleon inelastic cross section at y/s^N - 200 GeV is taken 
to (Tnn = 42 mb and the maximum nucleon density to po - 
0.17 nucleons/fm^. 



A. Partonic process for cc production 

Most studies of CNM effects on J /iff production [30] rely 
on the assumption that the cc pair is produced by a 2 — » 1 
partonic process where both initial particles are gluons carry- 
ing some intrinsic transverse momentum kr- The sum of the 
gluon intrinsic transverse momentum is transferred to the cc 
pair, thus to the J/tf> since the soft hadronisation process does 
not significantly modify the kinematics. This is supported by 
the picture of the Colour Evaporation Model (CEM) at LO 
(see [13] and references therein) or of the CO mechanism at 
[27]. In such approaches, the transverse momentum Pj of 
the J/ 1// is meant to come entirely from the intrinsic transverse 
momentum of the initial gluons. 

As just discussed, recent CSM-based studies [14, 15] have 
shown agreement with the PHENIX pp data [31] and the 
problematic Pj dependence of the LO CSM has been shown 
to be cured at Tevatron energies when going at or^ and [18- 
22]. Furthermore, e^e" J/t// + X^on cc NLO computations 
in the CSM [23, 24] leave now too small a room for a CO 
component which would support a 2 — > 1 production kinema- 
tics [27]. 

Parallel to this, intrinsic transverse momentum is not suf- 
ficient to describe the Pj spectrum of quarkonia produced in 
hadron collisions [13]. For Pj above approximately 2-3 GeV, 
one expects that most of the transverse momentum should 
have an extrinsic origin, i.e. the y/i/r's Pj would be balanced 
by the emission of a recoiling particle in the final state. The 
Jjijj would then be produced by gluon fusion in a 2 — > 2 pro- 
cess with emission of a hard final-state gluon. This emission 
has a definite influence on the kinematics of the 7/(/r produc- 
tion. Indeed, for a given Jjili momentum (thus for fixed y and 
Pt), the process g+g — > Jj^+g will proceed on average from 
gluons with larger Bjorken-x than § -H g — > cc — > Jjil/i+X). 
Therefore, they will be affected by different shadowing cor- 
rections. From now on, we will refer to the latter 2 — > 1 
scenario as the intrinsic scheme, and to the former 2 — » 2 as 
the extrinsic scheme. 

In the intrinsic scheme, we use the fits to the rapidity y 
and Pj spectra measured by PHENIX [31] in pp collisions at 
yfs^ - 200 GeV as inputs of the Monte-Carlo. The measure- 
ment of the Jlijj momentum completely fixes the longitudinal 
momentum fraction carried by the initial partons: 



particle introduces further degrees of freedom in the kinema- 
tics allowing more than one value of {x\,X2) for a given set 
{y, Pt) - which are ultimately the measured quantities as op- 
posed to {xi,X2) ■ The quadri-momentum conservation expli- 
citly results in a more complex expression of .51:2 as a function 
of (xi,y, Pt): 
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x-i nij y/sffffC y - 



(2) 



Equivalently, a similar expression can be written for 
xi as a function of {x2,y,PT)- In this case, mo- 
dels are mandatory to compute the proper weighting of 
each kinematically allowed {xi,X2). This weight is sim- 
ply the differential cross section at the partonic level 
times the gluon Parton Distribution Functions (PDFs), i.e. 
g(xi,fif)g(x2,fif)dcrgg^j/^+g/dydPTdxidx2. In the present 
implementation of our code, we are able to use the partonic 
differential cross section computed from any theoretical ap- 
proach. For now, we use the one from [14] which takes into 
account the i-channel cut contributions [32] to the basic CSM 
and satisfactorily describes the data down to very low Pt ■ A 
study using other matrix elements (LO CSM, NLO CEM,. . . ) 
is planned for future works. 
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FIG. 1: (Color online) comparison of the gluon shadowing 
parametiisations EKS98 [33], EPS08 [34], nDS [35] at LO and 
EPS09 [36] at LO in a Lead nucleus at Q- = 1.69 GeV^. Adapted 
from [36]. 



Xl,2 = 



-— exp i+y) = x'l 2(y, Pt), 



(1) 



with the transverse mass nij = -^M^ -1- P^, M being the JJij/ 
mass. 

On the other hand, in the extrinsic scheme, information 
from the data alone - the y and Pj spectra - is not sufficient 
to determine x\ and X2. Indeed, the presence of a final-state 



B. Shadowing 

To get the J I if/ yield in pA and AA collisions, a shadowing- 
correction factor has to be applied to the Jjij/ yield obtained 
from the simple superposition of the equivalent number of pp 
collisions. This shadowing factor can be expressed in terms 
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FIG. 2: (Color online) nuclear modification factor in dAu collisions, Rjau, at -^/Jnn = 200 GeV versus y in the intrinsic (top row) and 
extrinsic (lower row) schemes using the following gluon shadowing parametrisations: EKS98 [33] (left), EPS08 [34] (middle), nDSg [35] 
(right). The vertical lines indicate qualitatively the antishadowing peak and its shifts towards larger y in the extrinsic scheme. 



of the ratios Rf of the nuclear Parton Distribution Functions 
(nPDF) in a nucleon of a nucleus A to the PDF in the free 
nucleon. 

These parametrisations provide the nuclear ratios Rf of the 
PDFs: 

Rf(x, Q^) = , \ , i = q,q,8- (3) 

The numerical parametrisation of Rf(x, Q^) is given for all 
parton flavours. Here, we restrain our study to gluons since, at 
high energy, J I if/ is essentially produced through gluon fusion 
[13]. 

We shall consider three different shadowing parametrisa- 
tions for comparison: EKS98 [33], EPS08 [34] and nDSg [35] 
at LO. Recently, a new parametrisation, EPS09 [36], has been 
made available. It offers the possibility of properly taking 
into account the eiTors arising from the fit procedure. Yet, 
in the case of gluon nPDF, as illustrated by Fig. 1, the region 
spanned by this parametrisation is approximately bounded by 
both the nDS and EPS08 ones. However, we shall not consider 
the nDS parametrisation which shows such small shadowing 
coiTections that no significant yield coiTections are expected. 
We shall prefer to use nDSg as done in other works since it 
provides a lower bound of EPS09 in the antishadowing re- 
gion. Furthermore, the central curve of EPS09 is very close to 
EKS98. In this context, we have found it reasonable to limit 
our studies to EKS98, EPS08 and nDSg and to postpone to 



a further study the analysis of the error correlations and their 
impacts of y/f/r production studies'. 

Note that the spatial dependence of the shadowing has been 



Note that, as shown in Fig. 1 , the EPS08 parametrisation shows a strong 
gluon shadowing at very small .v due to the inclusion of forward-rapidity 
BRAHMS data in the tit. An explanation for such an effect has been based 
on the idea that, in this kinematic region that coiresponds to the beam frag- 
mentation region at large Feynman xp, one can reach the smallest values of 
the momentum fraction variable X2 in nuclei. It makes it possible to access 
the strongest coherence effects such as those associated with shadowing 
or alternatively the Color Glass Condensate (CGC). Nevertheless, as it has 
been argued in [37], although at forward rapidities one accesses the small- 
est X2 in the nuclear target, one simultaneously gets into the region of large 
xi of the projectile nucleon where energy conservation becomes an issue. 
Indeed, as stated in [38], at large xp (or xi ) one expects a suppression from 
Sudakov form factors, giving the probability that no particle be produced 
as Xf — > 1 , as demanded by energy conservation. In a pA collision, the 
multiple interactions of the nucleon remnants with the nucleus makes this 
less likely to occur and the suppression is expected to be stronger. In fact, 
factorisaton breaks down and the effective parton distributions in the pro- 
jectile nucleon then become con'elated with the nucleus target, see Eq. (15) 
of [38]. This effect should not be confused with gluon shadowing or other 
manifestations of coherence. Because of this, the strength of EPS08 gluon 
shadowing may be overestimated, since it includes data at large xp which 
were analysed assuming that the suppression was attributed only to a re- 
duction of the gluon distribution in the nucleus. However, we would like 
to emphasize that, while the gluon shadowing in EPSOS and EKS98 differ 
by a factor 3 at 2^ = 1.69 GeV", the difference is already less than 20% at 
=: M^j. which is the relevant scale for our analysis. 
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included in our approach, assuming an inhomogeneous sha- 
dowing proportional to the local density [39, 40]. 



C. The nuclear absorption 

The second CNM effect that we have taken into account 
concerns the nuclear absorption. In the framework of the 
probabilistic Glauber model, this effect refers to the proba- 
bility for the pre-resonant cc pair to survive the propagation 
through the nuclear medium and is usually parametrised by 
an effective absorption cross section cr^bs. Our results will be 
first shown for different values of cr^bs using the three afore- 
mentioned shadowing parametrisations. Afterwards, we will 
extract the values that provides the best fit to the PHENIX 
data. We note here that this effective cross section may also 
account for initial state effects, such as parton energy loss in 
the nuclear target. 



III. RESULTS FOR dAu COLLISIONS 

The J /iff suppression is usually characterised by the nuclear 
modification factor Rab, i-e., the ratio obtained by normalising 
the J /iff yield in ion collisions to the J /iff yield in proton colli- 
sions at the same energy times the average number of binary 
inelastic nucleon-nucleon collisions Ncii- 



dN- 



J/<l' 



AB 



(4) 



Any nuclear effect affecting J/iff production leads to a devia- 
tion of Rab from unity. 



A. RriAu vs rapidity: distribution shift 

In Fig. 2, we show R^au vs y according to both the ex- 
trinsic and intrinsic schemes. The results are displayed for 
four values of crabs for each of the aforementioned shadowing 
parametrisations and are compared with the PHENIX data [3]. 
The best fit result as performed in section V is also shown. 

The comparison between the three plots on the upper row 
and their corresponding plots on the lower one shows a strik- 
ing - but expectable - feature : the rapidity distribution in 
the extrinsic scheme is systematically shifted toward larger 
y compared to the intrinsic case. This is particularly visible 
when one focuses on the anti-shadowing peak, which we have 
indicated qualitatively with vertical lines. 

Such a shift is, in fact, not surprising at all. It simply re- 
flects the larger value of the gluon momentum fraction in the 
nucleus, X2, needed to produce a. J /iff when the momentum 
of the final state gluon is explicitly taken into account in the 
computations. 

As mentioned above, recent theoretical studies of J/^ pro- 
duction in vacuum, i.e. in pp collisions, support at low and 
mid Pj a partonic production mechanism as given by LO 
pQCD, namely gg — > J/iffg, as opposed to a 2 — » 1 process. In 



this context, we claim that this rapidity shift - evident for any 
shadowing parametrisation - is a feature of J/^ production 
in t/Au that should be systematically accounted for. Along 
the same lines of arguments, we shall focus in the following 
discussions on the results obtained in the extrinsic scheme 
(2^2 case), except for the extraction of crabs using the Rjau 
and Rcp results. 



B- RdAu vs centrality 

In Figs. 3 we present our results for R^au versus centrality, 
expressed as the number of collisions. We have taken into 
account the three shadowing parametrisations and four val- 
ues of (Tabs. One observes the effect of the impact parameter 
dependence of the shadowing -increasing for inner produc- 
tion points- which induces a progressive increase (resp. de- 
crease) vs A^coii in the backward (resp. forward) due to the 
anti-shadowing (resp. shadowing) effect. Indeed, for colli- 
sions with larger A^coib the J /if/ creation points are on average 
closer to the gold nucleus center where the shadowing is ex- 
pected to be stronger. For the same reason, the absorption 
suppresses the yield more strongly for larger A^coii- 

The overall effect (see Figs. 3) matches the trend of the 
PHENIX data [3], showing a decrease vs A^coU stronger in the 
forward region than in the backward. 

We also note that STAR collaboration has recently released 
a preliminary measurement of R^au in the region |y| < 0.5 
using the most central collisions (0 - 20%) [10]: RdAu - 
1.45 + 0.60. Higher statistics are however needed to draw 
conclusions from those data. 



C. RjAu vs transverse momentum 

It is important to note that, in order to predict the trans- 
verse momentum dependence of the shadowing corrections, 
one needs to resort to a model which contains an explicit de- 
pendence on Pj. Studies were earlier carried on using the 
CEM at NLO in [42]. However, due to the complexity inher- 
ent to the NLO code used, it was not possible to implement the 
impact parameter dependence of the shadowing, needed to re- 
produce, for instance, the centrality dependence [39, 40] as 
just discussed. Thanks to the versatility of our Glauber code, 
we can carry on such computation including such an impact 
parameter dependence as well as involved production mecha- 
nisms containing a non-trivial dependence on Pj- 

In Figs. 4, we show our results on R^au versus the transverse 
momentum. We emphasise that the growth of R^au is not re- 
lated to any Cronin effect, since it is not taken into account 
here. It simply comes from the increase of X2 for increasing 
Pj as given by Eq. (2). Hence, in the mid and forward ra- 
pidity region, xi goes through the anti-shadowing region and 
one observes an enhancement in RdAu- In the backward re- 
gion, where xo sits in anti-shadowing region, one only sees a 
decrease. A similar behaviour is obviously expected in Raa vs 
Pt as wUl be discussed in section IV C. 
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FIG. 3: (Color online) J/if/ nuclear modification factor in dAu collisions, RdAu, at -s/snn = 200 GeV versus the number of collisions for three 
rapidity ranges and for four values of the nuclear absorption (from top to bottom: cr„j^ =0, 2, 4 and 6 mb) using a) EKS98, b) EPS08 and c) 
nDSg in the extrinsic scheme. 
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FIG. 4: (Color online) J/if/ nuclear modification factor in dAu collisions, Rmu, at ■\fsNN = 200 GeV versus the transverse momentum for three 
rapidity ranges and for four values of the nuclear absorption (from top to bottom: Cabs =0, 2, 4 and 6 mb) using a) EKS98, b) EPS08 and c) 
nDSg in the extrinsic scheme. 
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FIG. 5: (Color online) J/i// central to peripheral nuclear modification factor in rfAu collisions, Rcp, at s/s^ = 200 GeV versus rapidity for 
three centrality ranges and for four values of the nuclear absorption (from top to bottom: a-„i,s =0, 2, 4 and 6 mb) using : a) EKS98, b) EPS08, 
c) nDSg. 



D. Rcp 

New results for the J/iJ/ from the 2008 dAu run with ap- 
proximately thirty times larger integrated luminosity than the 
2003 dAn results are emerging, with the first preliminary re- 
sult in terms of Rcp [11], 
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Those recent data show a significant dependence on the ra- 
pidity, with a visible suppression for the most forward points 
in the three centrality ranges (0-20%, 20-40% and 40-60%). 
In the negative rapidity region, which is expected to be domi- 
nated by large X2 contributions^, the data (see Fig. 5) show 
approximately no nuclear effects within the uncertainties, but 
a visible suppression is observed at mid and forward rapidi- 
ties. We also note that the suppression is stronger in most 
central collisions. As regards our results. Fig. 5 shows the ra- 
pidity dependence of the modification factor Rcp for the three 
centrality ranges using the EKS98, EPS08 and nDSg shado- 
wing parametrisations. Different break-up cross sections have 
also been considered, from to 6 mb. As announced, we have 
chosen to focus on the extrinsic scheme. Like for RdAu, this 
induces differences on the shadowing impact compared to an 



intrinsic scheme used in [28]. One observes on Fig. 5 that 
the trend of the data is reasonably described with a cr^bs in 
the range of 2-4 mb, while the most forward points seem to 
decrease more than our evaluation. However, the uncertain- 
ties affecting the present preliminary measurements precludes 
drawing firmer conclusions. We shall analysis this in more 
detail in section V. 



IV. RESULTS FOR CuCu AND AuAu COLLISIONS 
A. Centrality dependence^ 

In Fig. 6 and Fig. 7, we present the centrality dependence 
of the nuclear modification factor RauAu and RcuCu in the for- 
ward and mid rapidity regions. This has been computed for 
the three shadowing parametrisations and for five cr^t,^. 

As we have already observed in [9], Raa is systematically 
smaller in the forward region that in the mid rapidity region 
(see also the next section). The difference increases for more 
central collisions since we have used an impact parameter de- 
pendent shadowing. While this difference (approximatively 
20% for rather central collisions) matches well the one of the 
data when iXabs - 0, one would have to choose a large cr.^\,s if 
one wanted to reproduce the normalisation of the AuAu data, 
disregarding any effects of Hot Nuclear Matter (HNM). How- 



^ One recalls here that in the extrinsic scheme there is no one-to-one mapping 
between the rapidity and the momentum fraction of the initial gluon. Yet 
they are correlated. 



As announced, we focus on the extrinsic case in all the following discus- 
sions. 
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FIG. 6: (Color online) centrality dependence of (solid lines) and -/?™,^„ (dashed lines) for five values of the nuclear absorption (from 

top to bottom: (Tabs = 0, 2, 4, 6 and 8 mb) using a) EKS98, b) EPS08, c) nDSg, compared with the corresponding PHENIX data [2] 



O 

3 

o 



0.5 



extrinsic-p.|. EKS 




"CuCu 

Cu+Cu |y| < 0.35 
Cu+Cu |yN[1.2 2.2] 
syst = ± 8 % 

^ ^ g|obal I I I I 



20 40 60 80 100 
'^part 

(a) EKS98 



3 

o 
o 



0.5 



extrinsic-p.|. EPS 




"CuCu 

Cu+Cu |y| < 0.35 
Cu+Cu lyNJI.2,2.2] 
syst = ± 8 % 

^ ^ g|obal I I I I 



20 40 60 80 100 
Npar, 

(b) EPS98 



3 

o 

3 

o 



=^2r 0.5 



extrinsic-p.|. nDSg 




Cu+Cu |y| < 0.35 



Cu+Cu |yN[1.2,2.2] 
syst = ± 8 % 



20 40 60 80 100 
'^part 

(c) nDSg 



FIG. 7: (Color online) centrality dependence of J'^^aCiT'' (solid lines) and /?™,c„ (dashed lines) for five values of the nuclear absorption (from 
top to bottom: cr„i,j =0, 2, 4, 6 and 8 mb) using a) EKS98, b) EPS08, c) nDSg, compared with the corresponding PHENIX data [2]. 
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nish. Note that for a crabs in the range of 2-4 mb, a difference 
remains. 

While most of the above discussion is similar for CuCu col- 
lisions using the EKS and EPS shadowing, one observes a pe- 
culiar feature for nDSg. Indeed, one does not observe any 
effect. This comes for the very weak shadowing encoded in 
nDSg for Cu. This is illustrated on Fig. 8, where one sees that 
nDSg shadowing in Cu nucleus ends up to be very small. In- 
deed, while there is a moderate difference between the EKS 
shadowing in Cu and in Au, there is a significant difference 
for nDSg. 



FIG. 8: (Color online) comparison of the gluon shadowing 
parametrisations EKS98 [33] and nDSg [35] in lead and copper nu- 
clei at = 14 GeV^". 

"The plot has been generated by the nPDF generator 
http : //lappweb . in2p3 . f r/lapth/npd£generator 



ever, for such large cr^bs, surviving Jji^ from inner production 
points would be so rare that the difference between shado- 
wing effects at mid and forward rapidities would nearly va- 



B. Rapidity dependence 

We now discuss the rapidity dependence of the nuclear 
modification factor in the case of CuCu and AuAu colli- 
sions. It is displayed on Fig. 9 and Fig. 10 for the three 
shadowing parametrisations and for five (Jaha- As in [9], Raa 
slightly peaks at y = reducing the need for recombination ef- 
fects [43] which are particularly concentrated in the mid rapi- 
dity region and which could elegantly explain the differences 
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FIG. 9: Raa vs y for AuAu collisions for five cr.^f,s (from top to bottom: crnj, =0, 2, 4, 6, 8 mb) using a) EKS98, b) EPS08, c) nDSg in 4 
centrality bins. 




-2-1012-2-1012 -2-1012-2-1012 -2-1012-2-1012 

y y y y y y 

(a) EKS98 (b) EPS08 (c) nDSg 



FIG. 10: Raa vs y for CuCu collisions for five o-^bs (from top to bottom: a^i,, =0, 2, 4, 6, 8 mb) using a) EKS98, b) EPS08, c) nDSg in 4 
centrality bins. 



of Raa between the forward and mid rapidity regions. 

As we have noted in the previous section, shadowing effects 
exhibit naturally such a rapidity dependence. This happens for 
the three shadowing parametrisations we have used for AuAu, 
confirming that this is a feature, rather than an accident. This 
effect is, however, reduced when an absorption cross section 
is accounted for. It is widely accepted that HNM effects are 
responsible for an extra suppression. If this suppression - by 
the creation of a QGP for instance - is not coiTelated to the 
path of the cc on its way out of the nuclei, it may not damp 
down the difference between n-^™'"''^ and R<=^"'"\ as does a 
larger a-^b^. 

C. Transverse-momentum dependence 

We now move on to the discussion of the transverse mo- 
mentum dependence in the mid rapidity region, analysed both 
by the PHENIX [45] and STAR [44] collaborations. As an- 



nounced during the discussion of the dAu results, Raa versus 
Pt increases with Pj (see Fig. 11). In fact, the growth par- 
tially matches the trend of the PHENIX and STAR data. We 
should, however, mention here that there is no consensus for 
now on whether one should expect a nuclear modification fac- 
tor larger than one for Pj around 8 GeV as seems to indicate 
the pubHshed STAR results [44]. 

In the case of a confirmation of a non-suppression of 
at large Pj, one could say that it does not behave as the other 
hadrons, which are significantly suppressed in central Cu+Cu 
collisions and for increasing Pj (see [46] for tt" and [47] for 
"heavy-flavour" muons) In fact, the y/i/r seems to adopt a be- 
haviour closer to the one of prompt photons [49] than to the 
one of other (heavy-flavoured) hadrons. We also note a non- 
suppression of (f) meson in central Cu+Cu collisions [48]. 

If one goes further, one may want to extract information 
about the production mechanism at work here. Indeed, al- 
though the energy loss of a colored object in CNM is lim- 
ited to be constant, rather than scaling with energy, by the 
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FIG. 11: (Color online) J/ifi nuclear modification factor in CuCu collisions at -^Jivjv = 200 GeV versus Pt for the mid rapidity region for five 
values of the nuclear absorption (from top to bottom: cr^j, =0, 2, 4, 6, 8 mb) using : a) EKS98, b) EPS08, c) nDSg. Those are compared with 
the two STAR points [44] and the PHENIX preliminary results [45]. 



Landau-Pomeranchuk-lVIigdal effect [51], its magnitude per 
unit of length will be significantly larger for a CO than for a 
CS state propagating in the nuclear matter. This is especially 
relevant since, in the mid rapidity region, the cc pair hadro- 
nises outside the nucleus. This would naturally lead us to the 
conclusion that it is rather a colourless state than a coloured 
one which propagates in the nuclear matter. 





2 

""libs Xmin 


EKS98 Int. 
EPS08 Int. 
nDSg Int. 


3.2 ±2.4 0.9 
2.1+2^ 1.1 
2.2+2 3 13 


EKS98 Ext. 
EPS08 Ext. 
nDSg Ext. 


s.ril 1.1 

3.6+1^ 0.5 
3.0+1^ 1.2 



V. EXTRACTION OF BREAK-UP CROSS SECTION BY 
FITS ON rfAU DATA 

By comparing our results in the intiinsic and extiinsic ap- 
proaches, we have learnt that one of the consequences of this 
kinematical change implies a shift of the rapidity distribution. 
The latter is shifted as a whole to larger values of rapidity in 
the extrinsic case. As usual a break-up cross section, cr^bs, 
has to be accounted for to describe the normalisation of RdAu- 
In practice, it is fit from the data and then used in the descri- 
ption of nucleus-nucleus collisions. In view of the differences 
of the shadowing impact induced by one or the other kinema- 
tics, it is natural to wonder what the coiTesponding variations 
of (Tabs fit to the data are. 



A. Fitting Rm,, data 

For this purpose, we have used PHENIX measurements of 
^rfAu [3] in order to obtain the best fit of crabs for each of 
the shadowing parametrisation considered in both the intrinsic 
and extrinsic schemes. Based on the method used by PHENIX 
in [3] and [41], we have computed the;^'^ in the different cases, 
including both statistical and systematic errors. 

By using the data on RdAu versus rapidity, we have obtained 
the values of (Tabs given in Table I for each of the shadowing 
parametrisations and for both production schemes. The result- 
ing curves are shown on Fig. 2 for dAu, on Fig. 12 and Fig. 13 
for AuAu and CuCu. 

In general, fitting the rapidity with a constant (Tabs leads 



TABLE I: 1X21,5 extracted from fit of i?^^„ (All cross section in units 
of mb) 

to more than acceptable x^- The largest ones are obtained 
for the nDSg parametrisation; this confirms the impression 
that parametrisation with "significant" shadowing and/or anti- 
shadowing (EKS98, EPS08) are preferred by the data (see 
also Fig. 2). The only systematic one can really see is that 
the 3 (Tabs values extracted using the extrinsic scheme are 
larger than the coiTesponding ones obtained with the intrinsic 
scheme, as expected due to the increase of X2 in the extrinsic 
case compared to the intrinsic one. For the time being, the 
data are not precise enough to draw further conclusions from 
such fits. 



B. Fitting Rcp data 

We have also fit the new Rcp data with a constant (Tats in 
each rapidity regions. Our results are given in Table 11 along 
with the one of [28] based on a 2 — > 1 kinematics and using 
the EKS98 shadowing. We need to make clear at this point 
that the intrinsic scheme used in [28] slightly differs from the 
one we have used for instance to fit the dAu data as shown in 
the previous section. The difference appears at the level of the 
running of the scale of the shadowed gluon distribution and 
the invariant mass of produced system. However, this is not 
expected to modify the following discussion (see [8]). 

As can be seen in Table II, it appears that the strong sup- 
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FIG. 12: (Color online) R/y^ vs y for AuAu collisions for the best cr^bs as obtained from the fit to the dAu data using a) EKS98, b) EPS08, c) 
nDSg in 4 centrality bins. 
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FIG. 13: (Color online) i?^^ vs y for CuCu collisions for the best cr^bs as obtained from the fit to the dAu data using a) EKS98, b) EPS08, c) 
nDSg in 4 centrality bins. 



y < y = J > all >' 



EKS98 Int. [28] 


5.2 H 
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t 1.3 
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t 1.4 


N/A 


EKS98 Ext. 


2.5 H 


:0.5 
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tO.5 


4.8 = 


tO.7 


3.2 


±0.4 


EPS08 Ext. 


3.2 H 


= 0.5 


2.5 = 


tO.5 


3.1 = 


tO.6 


2.9 


±0.3 


nDSg Ext. 


1.4 H 
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tO.5 
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tO.7 
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TABLE II: cr;,bs and their errors" extracted from fit of the Rcp data 
averaged on three rapidity region as well as on the entire range. (All 
cross section in units of mb). 

"The errors quoted for the first line are extracted differently than ours. The 
errors shown here should only be compared for a given analysis. 



pression at forward rapidity and the lack of suppression at 
backward rapidity cannot be described using a fixed breakup 
cross section with EKS98 in the intrinsic scheme [28]. One 
also seems to observe an increase of cTabs with rapidity in our 
analysis (extrinsic) for EKS98 and nDSg, but a constant be- 
haviour cannot be ruled out. The increase is in any case softer 



when the final-state gluon momentum is taken into account. 
It seems that forward rapidity is maybe the most interesting 
region for such investigations. Interestingly, in the case of 
EPS08, the value we have extracted for the forward region is 
equal to the one for the backward region. 

As we have argued earlier, EPS08 can be used as a good 
indicator of the strongest shadowing reachable within the un- 
certainty of EPS09. From this viewpoint, the recent update 
of the intrinsic analysis, where [29] no increase of (Tabs is ob- 
served for the strongest shadowing of EPS09, confirms our 
findings. 

To get a more precise view on the situation, it is useful to 
plot the effective absorptive cross section as function of the 
rapidity, without averaging on the three rapidity regions. The 
result can be seen on Figs. 14 with our result in the extrin- 
sic scheme (red closed circles) and the one of [28] (blue open 
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(a) EKS98 (b) EPS08 (c) nDSg 



FIG. 14: (Color online) cr^i, versus y obtained by fitting the Rcp data for dAu, in the extrinsic scheme (in red closed circles) compared to the 
intrinsic scheme [28] (in blue open circles) using a) EKS98, b) EPS08 and c) nDSg. 



circles)"* in the intrinsic scheme. We have obtained an increas- 
ing absorptive cross section with rapidity, but the increase is 
much less pronounced than in the intrinsic case. In fact, in the 
EPS08 case, the increase does not appear statistically signifi- 
cant. 



VI. CONCLUSION 

Taking advantage of the probabilistic Glauber Monte-Carlo 
framework, JIN, discussed in [8, 9], we have (i) considered 
three different gluon shadowing parametrisations - EKS98, 
EPS08 and nDSg - taking into account a dependence on the 
impact parameter b and the momentum of the gluon recoil- 
ing against the J /if/, (ii) shown that the rapidity dependence of 
RdAu is shifted towards larger rapidities irrespective of the sha- 
dowing parametrisation, and particularly the anti-shadowing 
peak, (iii) shown that the anti-shadowing peak is reflected in a 
rise of the nuclear modification factor for increasing Pj, (iv) 
compared our results with the experimental measurements of 



the nuclear modification factors Rjau and Rcp from dAu col- 
lisions presently available at RHIC and extracted the favoured 
values of the cc absorption cross section in the nuclear mat- 
ter, and finally (v) shown that the effective absorption cross 
section increase at forward rapidity, obtained from the recent 
analysis of PHENIX Rcp data [28] in which the final-state 
gluon momentum is neglected (intiinsic case), is less marked 
when it is taken into account; that is in the extrinsic case. 
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